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Abstract 

We consider the tunnelling of charged spin-i fermions from a Kerr-Newman black hole and demon- 
strate that the expected Hawking temperature is recovered. We discuss certain technical subtleties related 
to the obtention of this result. 
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1 Introduction 

Semi-classical methods of modeling Hawking radiation as a tunnelling effect were developed over the past 
decade and have garnered a lot of interest pQ-[55]. The earliest work with black hole tunnelling was done by 
Kraus and Wilczek [TJ, and was refined by various researchers [31 H] - From this approach an alternative way 
of understanding black hole radiation emerged, fn particular one can calculate the Hawking temperature in a 
manner independent of traditional Wick Rotation methods or the original method of modelling gravitational 
collapse |29j . Tunnelling provides not only a useful verification of thermodynamic properties of black 
holes but also an alternate conceptual means for understanding the underlying physical process of black 
hole radiation. For scalar field emission it has been shown to be very robust, having been successfully 
applied to a wide variety of interesting and exotic spacetimes, including the Kerr and Kerr-Newman cases 
[TOl [nj [TJ], black rings [TJ], the 3-dimensional BTZ black hole [13], the Vaidya spacetime [18], other 
dynamical black holes [12], Taub-NUT spacetimes pH] , Godel spacetimes |22j . and dynamical horizons [19 . 
Tunnelling methods have even been applied to horizons that are not black hole horizons including those with 
cosmological horizons [5] .[5] .[2"5] .and Rindler Spacetimes [3] , [TJ] , (55] for which it has been shown the Unruh 
temperature [30) is in fact recovered. 

In general tunnelling methods involve calculating the imaginary part of the action for the (classically 
forbidden) process of s-wave emission across the horizon, which in turn is related to the Boltzmann factor 
for emission at the Hawking temperature. Different approaches exist for calculating the imaginary part of 
the action for the emitted particle. The first black hole tunnelling method developed was the Null Geodesic 
Method used by Parikh and Wilczek [3] , which followed from the work of Kraus and Wilczek [T] . The other 
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approach to black hole tunnelling is the Hamilton- Jacobi Ansatz used by Agheben et al [7], which is an 
extension of the complex path analysis of Padmanabhan et al [4] . Recently we have extended black hole 
tunnelling to include the emission of spin 1/2 fermions [21], demonstrating that fermion tunnelling from 
both a Rindler horizon and a generic non-rotating black hole recovers the expected results for temperature. 
While perhaps not surprising, the result is non-trivial insofar as fermionic vacua are distinct from bosonic 
vacua, and can lead to distinct physical results [3T] . This work has been extended to describing fermionic 
tunnelling across horizons in other spacetimes such as the BTZ black hole [26] . and dynamical horizons |27j . 

All tunnelling approaches use the fact that the WKB approximation of the tunnelling probability for the 
classically forbidden trajectory from inside to outside the horizon is given by: 

r oc exp(-2Im/) (1) 

where / is the classical action of the trajectory to leading order in H (here set equal to unity). Where these 
methods differ is in how the action is calculated. For the Null Geodesic method the only part of the action 
that contributes an imaginary term is J r ° ut p r dr, where p r is the momentum of the emitted null s-wave. Then 
by using Hamilton's equation and knowledge of the null geodesies it is possible to calculate the imaginary 
part of the action. For the Hamilton-Jacobi ansatz it is assumed that the action of the emitted (scalar) 
particle satisfies the relativistic Hamilton-Jacobi equation. From the symmetries of the metric one picks an 
appropriate ansatz for the form of the action and plugs it into the Relativistic Hamilton-Jacobi Equation 
to solve. (For a detailed comparison of the Hamilton-Jacobi Ansatz and Null-Geodesic methods see [T4j). 
The Hamilton-Jacobi Ansatz came from applying the WKB approximation to the Klein-Gordon equation. 
To lowest order in WKB this results in the Hamilton-Jacobi equation. For tunnelling of spin 1/2 particles 
it can be shown that applying the WKB approximation to the Dirac equation instead of the Klein-Gordon 
equation yields the tunnelling probability for fermions [21]. This was the first time that the tunnelling 
approach had been used to model spin 1/2 particles. 

In this paper we extend the tunnelling method to model charged spin 1/2 particle emission from rotating 
black holes. To this end we apply the fermion tunnelling method to the Kerr-Newman black hole for both 
massless and massive charged particle emission. This extension introduces some non-trivial technical features 
associated with the choice of 7 matrices. We confirm that spin 1/2 fermions are emitted at the expected 
Hawking Temperature from rotating black holes, providing further evidence for the universality of black hole 
radiation. 
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2 Charged Spin 1/2 Particle Emission From Kerr-Newman Black 
Holes 



We will consider particle emission from the Kerr-Newman solution. The Kerr-Newman metric and vector 
potential are given by 

ds 2 = - f(r, 0)dt 2 + — — - 2H(r, 9)dtd(f> + K(r, 6)d(j) 2 + E(r, 6)d9 2 

g(r,o) 

Aa = -^P) a - a sin 2 0(#) a ] (2) 
E(r) 



A(r) - a 2 sin 2 ■ 
A(r) 



/M) 
<?M) 
ff(r,0) 



EM)' 

a sin 2 6»(r 2 + a 2 - A(r)) 
EM) 



E(r, 0) 
E(r,0) = r 2 +a 2 cos 2 
A(r) = r 2 + a 2 + e 2 - 2Mr 

Since the tunnelling method is not applicable to extremal black holes [14] , we will assume a non-extremal 
black hole so that M 2 > a 2 + e 2 . Consequently there are two horizons located at r± = M ± \J M 2 — a 2 — e 2 . 
It is convenient for our calculations to work with the function F(r, 9) — — (g 11 ) where 



nr, 0) = f(r, 9) + tL^L = _ . (3) 



H 2 (r,9) _ A(r)S(r,0) 
K(r, 9) ~ (r 2 + a 2 ) 2 - A(r)a 2 sin 2 9 

and where the angular velocity at the black hole horizon is 



H(r + ,6) a 

= W^ff) = (4) 

We will only show the calculation explicitly for the spin up case; the final result is also the same for 
the spin down case as can be easily shown using the methods described below. In the non-rotating case 
a statistical argument was used to justify the assumption that overall a zero angular momentum state is 
maintained for fermion emission, because as many particles with spin pointing radially outward (spin up) 
would be emitted as particles with spin pointed radially inward (spin down). This argument is still valid in 
the rotating case: the statistical distribution of spins in the fermion emission spectrum should not alter the 
angular momentum of the black hole. 

The Dirac equation with electric charge is: 



^(^-|^+^ = (5) 
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where: 



D„ 



a. 



d u + a 



(G) 
(7) 

(8) 



The 7 M matrices satisfy {7 A ',7 iy } = 2g^ v x 1. We choose a representation for them in the form: 



7 = 



7 2 + 



y/K^ff) \ y/F(r,0)K(r,9) J 
where the 7°'s are simply the following chiral 7's for Minkowski space 



(9) 



7° = 



7 2 = 



and the a' s are the Pauli Matrices 
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and we denote £,y^ for the eigenvectors of a 3 . Note that 
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7 = 27 7 7 7 



FKYj 



/FK 







'FK 



is the resulting 7 s matrix. 

The spin up (i.e. +ve r-direction) ansatz for the Dirac field, has the form: 



A(t,r,M)£ t 
B((t,r,M)£ T 

A(i,r,M) " 






exp 



-Mt,r,6,(j>) 



exp 



-I r (t,r,6,<l>) 



(10) 



(11) 



(12) 



(13) 



In order to apply the WKB approximation we insert the ansatz (|13p for spin up particles into the Dirac 
Equation. Dividing by the exponential term and multiplying by H the resulting equations to leading order 
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Note that although A, B are not constant, their derivatives - and the components f2 M - are all of order 0(h) 
and so can be neglected to lowest order in WKB. 

When ra^O equations (fT4]) and (fl"6|) couple whereas when to = they decouple. We employ the ansatz 



(18) 



and insert it into equations (|14M17[) (where we consider only the positive frequency contributions without 
loss of generality). To simplify the expressions we expand the equations near the horizon and find 
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In the massless case it is possible to pull , 1 out of equations (fTl?|) and (f2"Tj) . making these equations 

ys(r + ,e) 



independent of 9. Furthermore, equations (|20| and ([22]) have no explicit r dependence. From this we can 
conclude that near the black horizon it is possible to further separate the function W 

W{r,6) = W{r) +6(0) 

and we see that equations (f2"0"|) and (j2"2"|) both yield the same equation for regardless of A or B. 
Equations p^|) and (12111 then have two possible solutions 



A = and W'(r) = W' + {r) 
B = and W'(r) = W'_{r) 



(E-Q H J-^)(rl + a 2 ) 



A r (r+)(r - r+) 

_qer+ 

r'i+a 



- { E-n H J-^)(r\ 



&r{r+)(r - r+) 



where the prime denotes a derivative with respect to r and W + /_ corresponds to outgoing/incoming solutions. 
Since we are only concerned with calculating the semi-classical tunnelling probability, we will need to multiply 
the resulting wave equation by its complex conjugate. So the portion of the trajectory that starts outside 
the black hole and continues to the observer will not contribute to the final tunnelling probability and can be 
safely ignored (since it will be entirely real). Therefore, the only part of the wave equation that contributes 
to the tunnelling probability is the contour around the black hole horizon. For a visual representation of the 
deformation of the contour see Figure[T] This contour differs somewhat relative to conventions Padmanabhan 
defines |4J, in which contours in the upper half plane are selected for both ingoing and outgoing particles. 
For the emission integral he multiplies his equivalent of W + (dSo/dr) by a minus sign "where the minus sign 
in front of the integral corresponds to the initial condition that dSo/dr > at r = r% < ro" [4]. Instead we 
choose the mathematically equivalent convention that the outgoing contour is in the lower half plane and so 
do not multiply by a minus sign. 

The probabilities of crossing the horizon in each direction are proportional to 

Prob[out] oc exp[-2Imi] = exp[-2(ImM / + +Im9)] (23) 
ProbH cx exp[-2Im7] = exp[-2(ImW_ +Im6)] (24) 

To ensure that the probabilities are correctly normalized so that any incoming particles crossing the 
horizon have a 100% chance of entering the black hole we need to divide each equation by (|24p. From 
this the probability of going from outside to inside the horizon will be equal to 1 and this implies that the 
probability of a particle tunnelling from inside to outside the horizon is: 

^ Prob[out] exp[-2(ImVF+ +Im6)] . AT „, . 

r oc ^ ,\. ' = n ) T — — — = cxp -4ImVF + (25) 

Prob m exp -2(ImWl +Im6) 1 +J y ' 



Solving for W+ yields 

W+(r) 



(E-n H J-^)(rl+a 2 ) 
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Figure 1: Diagram of contours between black hole and observer for both outgoing and incoming trajectories 
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and after integrating around the pole (and dropping the + subscript) we obtain 

m {E-Sl H J-^)(r\ + a 2 ) 



W = 



' + 



2r i - 2M 



The resulting tunnelling probability is 



r = exp[-27r '+^1 (E - Q H J - -4^)} 



giving the expected Hawking temperature 



= 1 r. M 1 (M 2 -q 2 -e 2 )^ 
H 2vr r\ + a 2 2vr 2M(M + (Af 2 - a 2 - e 2 ) 3 ) - e 2 

for a charged rotating black hole. 

In the massive case equations (|19|l and (j2"Tj) no longer decouple and analysis of the tunnelling is more 
subtle. We begin by eliminating the function W r (r, 9) from these two equations so we can find an equation 
relating A and B in terms of known quantities. Subtracting Bx (|2"T]) from Ax (JTSJ) gives 

2AB(E-{l H J-^±_) 

= = t=? —+mA 2 -mB 2 = (28) 

y/F r (r + ,6)(r-r+) 

,A o „ <Z er + s/^- 



= mVF r (r+, - r+)(-) 2 + 2(E - U H J - -*_±_)(_) 



iy/F r (r + ,6)(r-r+) (29) 



and so 



A —{E flff J ^) ± ^{E n H J - ^) 2 + m*F r (r + ,9)(r r+) 
B m v / F r (r + ,9)(r - r+) 



where 



-(25 - fi* J - ^) ± /(£ - fi H J - + m^F r (r + , 9)(r - r+) \ , 

lim I V - + = <^ (30) 

r-r+ I m v /Fr(r + ,0)(r-r + ) 1 1 ^ 



-oo 



for the upper/lower sign respectively. 

Consequently at the horizon either 4 — > or — > — oo, i.e. either ^4 — > or B — > 0. For A — > at the 
horizon, we solve (|2ip in terms of m and insert into (|19p , obtaining 



v / F r (r + ,% r (r + ,0)(r - r+) (l - #) 



^(r, ) = „, , ==; — ( 31 ) 
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Note that the ^-dependence drops out of this expression, i.e. 




since -g is zero at the horizon the result of integrating around the pole is the same as in the massless case. 
For B — > we can simply rewrite the expression (l3~Tj) in terms of to get 




Again, since the extra contributions vanish at the horizon, the result of integrating around the pole for W 
in the massive case is the same as the massless case and we recover the Hawking temperature (|27p for the 
Kerr-Newman black hole. 

The spin-down case proceeds in a manner fully analogous to the spin-up case discussed above. Other 
than some changes of sign the equations are of the same form as the spin up case. For both the massive 
and massless cases the temperature (I27[) is obtained, implying that both spin up and spin down particles are 
emitted at the same rate. 

3 Technical Issues 

With rotating spacetimes the choice 7 matrices is quite relevant, not only for ease of calculation but also for 
tractability. In order to demonstrate this we will repeat the calculation for a different (and less convenient) 
choice of 7 matrices. We will also set the charge q of the emitted particles to zero for simplicity. 
Consider the choice 



where we use the same chiral 7° matrices (ITU)) as before. This choice satisfies the correct anti-commutation 
relations {7^, 7"} = 2g^ , and corresponds to a different choice of tetrad basis for the metric. 
Naively choosing the same ansatz as before 




(32) 



^ T (t,r,M) 
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exp 




we find upon insertion into the (chargeless) Dirac equation §5§ we obtain 
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(33) 
(34) 
(35) 
(36) 



-Et + J(j> + W ) we find from §3§ and 



Repeating the same kind of analysis as before (using the ansatz / 
(j35|) that 

w r± = ~- E 



Since / does not vanish at the horizon (except when sin# = 0), this expression does not have a simple 
pole at the horizon. It is not possible to solve the expression for arbitrary 6, and the calculation becomes 
intractable. This situation is analogous to what happens in the scalar field case if one naively applies the 
null geodesic method to a rotating black hole by trying to force <j> to be constant (id d(f) = 0), as previously 
demonstrated [14] . 

In order to understand this issue in more detail it is useful to examine the similarity transformation 
between 7^ and 7^. We find that: 



7" = sys -1 , 



when: 



s = 



al - ba 2 




for all [i 





al + ba 2 



where 



a = 



\ 




\ 




(37) 



The transformation S is similar to a Lorentz boost in the 4> direction. Applying it to the spin up ansatz 
used previously we find 



ip r (t,r,e,<t>) = 



Aa 
-Aib 
Ba 
Bib 



exp 



(38) 



As r — > 00 we see that a — > 1 and b — > 0, yielding the same spin up ansatz in this limit. Inserting (|38[) into 
the (chargeless) Dirac equation ([5]) and following the same procedure as before results in the same expression 
(|27p for the temperature. This is not surprising since all we have done is applied a similarity transformation 
to the Dirac equation, and we shall not repeat the (somewhat more tedious) calculations here. Our point 
is to emphasize the importance of choosing an appropriate ansatz for a given choice of 7 matrices. 
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4 Conclusions 



We have successfully extended our approach of fermion tunnelling to model the emission of charged fermions 
from a rotating charged black hole. The analysis yields the expected temperature (|2T[) . consistent with black 
hole universality. However there are subtle technical issues involved with choosing an appropriate ansatz 
for the Dirac field consistent with the choice of 7 matrices, and failure to make such a choice leads to a 
breakdown in the method. 

Computing corrections to the tunnelling probability by fully taking into account conservation of energy 
will yield corrections to the fermion emission temperature. In various scalar field cases this is inherent in the 
Parikh/Wilczek tunnelling method [3], [H]-[inj and can be incorporated into the Hamilton- Jacobi tunnelling 
approach [8] . Another avenue of research is to perform tunnelling calculations to higher order in WKB (in 
both the scalar field and fermionic cases) in order to calculate grey body effects. It is also worth investigating 
the possibility of calculating a density matrix for the emitted particles from a tunnelling approach in order 
to calculate correlations between particles. Work is in progress in these areas. 
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